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Abstract 

This report presents the results of a pretest structural analysis of an 
all-composite stitched/RFI (resin film infusion) wing stub box which was 
designed and fabricated by the McDonnell Douglas Aerospace Company. 
Geometrically nonlinear structural responses of the wing stub box were 
predicted by using the finite element analyses and a global/local approach 
in which the global model contains the entire test article while the local 
model contains a large nonlinearly deformed region in the upper cover of 
the wing stub box. The wing box test article includes the all-composite 
wing stub box, a metallic load-transition box and a metallic wing-tip 
extension box. The two metallic boxes were connected to the inboard and 
outboard ends of the composite wing stub box, respectively. The metallic 
load-transition box was attached to a steel and concrete vertical reaction 
structure. In the global analysis, an upward load was applied at the tip of 
the extension box to induce bending of the wing stub. Global analysis found 
that an upper cover region, which contains three stringer runouts, exhibits 
large nonlinear deformations. Hence, a local model refined in the 
nonlinearly deformed region was created to predict more accurate strain 
results near stringer runouts. Numerous global and local analysis results 
such as deformed shapes, displacements at selected locations, and strains at 
critical locations are included in this report. 

Introduction 

The purpose of this report is to document the pretest structural 
analysis results for an advanced stitched/RFI (resin film infusion) 
graphite/epoxy wing stub box. This advanced stitched/RFI graphite/epoxy 
wing stub box, representing the inboard portion of a civil-transport- 
aircraft high-aspect-ratio wing stub box, was designed and manufactured 
by McDonnell Douglas Aerospace (MDA) Company under the support of 
NASA’s Advanced Composites Technology (ACT) Program. The 
innovative stitched/RFI process used to fabricate this all-composite wing 
stub box has the potential for reducing manufacturing costs while 
producing damage tolerant composite aircraft primary structures. 



A preliminary global analysis was performed by McDonnell Douglas 
engineers[l]. More refined nonlinear analyses of the global model and a 
local model were performed by NASA Langley Research Center engineers 
to support the structural wing-stub-box tests at NASA Langley. The global 
model contains the entire test article which includes an all-composite stub 
box, an inboard metallic load transition box and an outboard metallic wing- 
tip extension box. The local model contains an upper cover region 
surrounding the location of three stringer runouts (where a stiffener 
terminates at a rib). It was found that this region exhibits large nonlinear 
deformations. Hence, a more refined mesh was used in the local model to 
obtain more accurate stress analysis results. Displacements obtained by the 
global model analysis were applied to the boundaries of the local model. 
Deformed shapes, displacements at selected locations, and strains in critical 
regions generated from the global and local analyses were used for 
instrumenting the wing stub box. Only the analysis results are presented in 
this report. These analytical results have been correlated with test results 
and the correlation will be presented in another paper [2]. 

Wing Stub Box Test Article 

A photograph of the wing stub box test article is shown in Figure 1 
and the dimensions of the wing stub box are provided in Figure 2. The 
composite wing stub box is about 12 ft long and 8 ft wide, and its 
maximum depth at the root is about 2.3 ft. The wing stub box test article 
includes the all-composite stub box, the inboard transition box (made of 
steel and aluminum), and the outboard extension box (made of steel). The 
all-composite stub box weighs approximately 1200 lbs while the entire 
wing stub box test article weighs about 7600 lbs. In its test configuration, 
the transition box was attached to a steel and concrete strongback (test wall) 
at the NASA Langley Structural Mechanics Test Laboratory. The extension 
box is a load introduction structure and the ultimate design load, a vertical 
load at the tip of the extension-box front spar, is 166,000 lbs. 

The upper-cover-panel construction of the all-composite wing stub 
box shown in Figure 3 contains a skin panel, ten blade stiffeners, five ribs, 
two metal angles, two spar caps, and has a 22.5 in. by 12.5 in. oval access 
door cutout. The upper cover skin, flanges, and blade stiffeners were 
constructed from AS4/3501-6 graphite/epoxy composite with repeating 
sublaminates (stacks). The layup sequence of a stack is [45/-45/0/90/0/- 
45/45] and the nominal total thickness of a stack after curing is 0.058 in. 
The upper skin thickness varies from five stacks (0.29 in.) to ten stacks 
(0.58 in.) as shown in Figure 4. All the blade stiffeners of the upper cover 
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have the same thickness and each has eight stacks (0.464 in.). The upper- 
cover panel was fabricated by first assembling and stitching all the dry 
preform components together and then using a resin film infusion process 
to infuse resin into the preform. Each dry preform component such as the 
skin or the stiffener was stitched individually before assembling the upper- 
cover-panel. The stiffener flanges are stitched to the skin so that no 
mechanical fasteners are required. This reduces the manufacturing costs. 
However, at the stiffener runout locations, fasteners were installed after the 
RFI process to prevent skin stiffener debonding at these locations. 
Moreover, the access door cutout was reinforced by a composite land ring 
which was bolted to the panel skin using 0.3125 inch diameter bolts. To 
further examine the region between Ribs 6 and 8, which contains three 
stringer runouts and two metal angles as shown in Figure 3, a local refined 
nonlinear analysis was carried out. 

The interior of the wing stub box is shown in Figure 5 in which the 
five ribs and the two spar webs are made of conventional AS4/3501-6 
prepreg material. The spar webs have a constant thickness of 0.31 in. and 
the rib webs have a constant thickness of 0.15 in. Moreover, the ribs and 
spars are stiffened with composite stiffeners to prevent buckling. The 
lower skin is made of stitched/RFI graphite/epoxy material with IM7 fibers 
in the 0-degree fiber direction and AS4 fibers in all the other fiber 
directions. The thickness of the lower skin, (which is thicker than the 
upper skin), ranges from .33 in. to .82 in. The high moduli of IM7 fibers 
and thicker skin result in smaller strains for the lower skin; therefore no 
failure is expected in the lower skin region and the results in this report 
relate mainly to the upper cover panel. 

Material Properties and Allowables 

The equivalent AS4/3501-6 material properties for the upper skin 
panel laminates used in the analyses are 

E x = 8.17 Msi 
E y = 4.46 Msi 
G_. = 2.35 Msi 
= 0.459 

and the equivalent AS4/IM7/3501-6 material properties for the lower skin 
panel laminates used in the analyses are 
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E x = 11.85 Msi 
E y = 4.55 Msi 
= 2.57 Msi 
= 0.409 

Note that the x-direction, which is parallel to the rear spar, is coincident 
with the 0-degree fiber direction and the y-direction is coincident with the 
90-degree fiber direction. The ultimate compression strain allowable in 
the x-direction of the undamaged upper skin laminate is 9330 and the 
0.3125 inch diameter filled hole B-base compression strain allowable is 

8100 jie [1]. 


Global and Local Finite Element Models 


The global finite element model shown in Figure 6 was used to 
determine the global responses of the complete test article and this global 
model contains 4408 quadrilateral elements (CQUAD4), 99 triangular 
elements (CTRIA3), 1308 beam elements (CBEAM), 798 rod elements 
(CONROD) and 742 rigid bar elements (RBAR). The finite element 
analyses were performed using MSC/NASTRAN [3], V68. The upper skin 
lower skin, rib webs and spar webs were modeled as plate elements 
(CQUAD4 and CTRIA3), and the stiffeners were modeled as beam 

elements (CBEAM). A total of 5,266 grid points were used in the global 
model. 

A finer finite element mesh was used around the access door cutout 
to better define the high stress concentration in that region. Grid points of 
each composite stiffener were positioned along its centroidal axis. 
Therefore, the stiffener elements, modeled as beam elements, were offset 
from the skin. These stiffener beam elements were rigidly connected to the 
skin by relatively stiff CBEAM elements. The finite element mesh of the 
wing stub box interior region is shown in Figure 7 in which all the rib 
webs including cutout holes were modeled. Both buckling analysis and 
geometrically nonlinear analysis of the global model were performed by 
using MSC/NASTRAN and the solution sequences used for buckling 
analysis and nonlinear static analysis were SOL 105 and SOL 106 
respectively. 

As mentioned earlier, the local finite element model contains an 
upper cover region surrounding three stringer runouts as shown in Figure 
8. The location of the local model on the stub box is clearly shown in 
Figure 9 in which the local model is superposed on the global model. 
Structural details near a typical stiffener runout region are shown in Figure 
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10. Note that the tapering of the blade stiffener flange and the blade web 
in the runout region and a small gap existing between the stiffener flange 
and rib flange were also modeled in the local model. The blade stiffener 
web is attached to the rib web, but it is not connected to the rib flange. The 
local finite element model shown in Figure 8 has 2642 quadrilateral shell 
elements, 8 triangular elements and 2738 nodes. Displacement boundary 
conditions of the local model were obtained through spline interpolation of 
the global displacement results. NASA Langley’s COMET (Computational 
Mechanics Testbed) finite element code [4] was used for the nonlinear 
analysis of this local model. 

The global model was originally created by engineers at the 
McDonnell Douglas Aerospace Company and some modifications were 
performed by engineers at NASA Langley. PATRAN [5] was used to 
modify the global model and was also used to create the local model. 
Results from both the global and local analyses were postprocessed using 
PATRAN. 


Analysis Results 


Global analysis results 

A buckling analysis of the wing-stub-box global model was 
performed first and the critical buckling load was found to be 11.2 % 
higher than the ultimate design load of 166,000 lbs with the buckling mode 
shape shown in Figure 1 1 . The second analysis performed was a 
geometrically nonlinear NASTRAN analysis and the predicted deformed 
shape of the whole test article at ultimate load is shown in Figure 12. The 
jack loading displacements predicted at ultimate load were u= -0.2169 in, 
v=0.1629 in, and w= 13.4887 in. The predicted relationship between the 
load and the vertical displacement at the loading point is approximately 
linear as shown in Figure 13. The close-up view of the deformed shape of 
the composite wing stub box is shown in Figure 14. Note that the wavy 
deformation of the upper skin is likely caused by a lack of longitudinal 
stiffener support near the access door cutout and its two adjacent outboard 
bays as shown in Figure 3. Figure 15 illustrates the out-of-plane (z- 
direction) deflections in the deformed region along two lines, labeled A-A 
and B-B, parallel to the major axis of the elliptic cutout. Line B-B, located 
at 15.875 in. from Line A-A, is in location which sees less wavy 
deformation and the difference of the out-of-plane displacements between 
these two lines provides a good measurement of the nonlinearity 
developing in the two bays outboard of the access door. 
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The predicted out-of-plane displacements at six LVDT (Linear 
Voltage Displacement Transducer) locations on the lower skin are plotted 
in Figure 16. Locations 1 to 3 are on the rear spar while Locations 4 to 6 
are on the front spar (see insert in Figure 16 and the Douglas Aircraft 
Company Tasks Assignment Drawings TAD Z7944503). All the out-of- 
plane displacements are seen to be linearly related to the applied load. 

All the strain results presented for the global model are based on the 
NASTRAN element coordinate system. The element x-axis is coincident 
with the global x-axis except in the access door cutout region where the 
element x-axis is tangent to the edge of the cutout (in the circumferential 
direction). Strain contour plots (e^) of the upper composite skin panel 
are plotted in Figures 17-19. High strains occur at the edge of the access 
door cutout. Close-up strain contour plots for the skin side, mid-surface 
and stiffener side of the upper cover in the cutout region are in Figures 20 
to 22. High strains are predicted in the circumferential direction of the 
cutout: 9100 fie at the skin side, 6650 fie at the stiffener side, and 6770 fie 
at the midsurface. 

Predicted strains at all strain gage locations for 100% and 50% of 
the ultimate load are listed in Tables I to IV in which the strain gage 
locations and the completed term of every abbreviation can be found in 
Douglas Aircraft Company Tasks Assignment Drawings TAD Z7944503. 
Note that the strain gage patterns for the upper and lower skins are shown 
in Figures 23 and 24. 

Strain versus load plots are generated for some strain gages located 
in high strain regions including these gages near the access door cutout 
(§ a gcs 78, 79, and 612) and in the highly nonlinear deformed regions 
(gages 63, 64, 67, and 68). Locations of these gages are shown in Figure 
25 and predicted strain versus load plots at these gage locations can be 
found in Figures 26 to 30. 

Local analysis results 

A geometrically nonlinear analysis was performed with the local 
model using the COMET code and the displacement boundary conditions 
obtained from the nonlinear global analyses. These displacement boundary 
conditions were applied on the four edges of the local model and the root 
of the middle rib. The predicted deformed shape of the local model at the 
ultimate loading condition is shown in Figure 31. Note that the skin is 
deformed to the stiffener side which may be due to the lack of longitudinal 
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stiffener support in one of the skin bays. The skin strain contour plots at 
the ultimate load, as derived from the nonlinear analysis, are plotted in 
Figures 32 to 34 at the skin side, stiffener side, and mid-surface. The 
strains in the skin side are much higher than in the stiffener side due to the 
significant skin bending deformation shown in Figure 31. 

Figures 35 to 37 display the variation of the strain in the x-direction 
with load for the three gaps between the stiffener flange and rib flange, 
(see Figure 10), located at the stiffener runouts, (labeled #1 to #3 in Figure 
8), for the skin side, mid-surface, and stiffener side locations. Predicted 
lateral strains (in y-direction) in the rib webs are plotted in Figures 38-40. 
Significant rib-web bending was predicted at the runouts where the 
stiffener webs attach to the ribs. 

Predicted strains at all strain gage locations in the local model region 
for 100% and 50% of the ultimate load are listed in Tables V and VI. Note 
that the strain gage locations can be found from Figure 23. 

The existence of the stiffener runouts and the large skin bending 
deformations may induce the stiffener-skin separation loads, defined as the 
vertical or normal stress resultants between the blade and the skin. The 
blade stiffeners in the local model are numbered from 1 to 4 as shown in 
Figure 8. The separation load of the front spar may not be an issue 
because the front spar cap is mechanically fastened to the skin, thus only 
the separation load results for blade stiffeners are presented and plotted in 
Figures 41 to 44. The horizontal axes of these plots are the distance 
measured from an edge of the local model in the positive x-direction (see 
insert of each plot). These results can be compared with blade stiffener 
separation allowables [6]. The ultimate failure load of blade stiffener 
separation tests performed in Reference 6 is about 2500 lbs/in which is 
much higher than the predicted separation loads here. Therefore, no 
stiffener-skin separation failure is expected. 

Concluding Remarks 

This report documents the global and local geometrically nonlinear 
finite element analysis results, generated by engineers at NASA Langley, of 
the McDonnell Douglas Stitched/RFI wing stub box. These analyses were 
performed by the Computational Structures Branch of NASA Langley 
Research Center. Results have been used for strength prediction of the stub 
box and determination of the instrumentation and gage locations of the stub 
box test article. When the actual structural tests are performed on the stub 
box, analysis results will be correlated with the test data. 
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Table I. Strain prediction from global model at ultimate load of 166,000 lbs 
(including all gages on external upper skin). 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

2 

0.0 

ext 


-3465.4 

196 

4 

0.0 

ext 

essshhi 

-3758.9 

194 

5 

0.0 

ext 

QisS^Si 

-5032.7 

65 

6 

45.0 

ext 


-1761.4 

65 

7 

90.0 

ext 


1821.7 

65 

11 

0.0 

ext 


-4978.8 

1359 

12 

0.0 

ext 

OJBSSSHB 

-4978.8 

1359 

15 

0.0 

ext 


-4556.0 

614 

16 

0.0 

ext 



614 

17 

0.0 

ext 

hkessm 


71 

18 

0.0 

ext 


-4842.3 

414 

20 

0.0 

ext 


-4953.4 

437 

22 

0.0 

ext 

O 

-5194.9 

371 

25 

0.0 

ext 


-4210.0 

101 

26 

0.0 

ext 

Upr.Skin 

-4433.7 

403 

27 

45.0 

ext 

Upr.Skin 

-805.2 

403 

28 

90.0 

ext 

Upr.Skin 

1395.3 

403 

32 

0.0 

ext 

Upr.Skin 

-4085.5 

401 

37 

0.0 

ext 

Upr.Skin 

-4750.7 

474 

39 

0.0 

ext 

Upr.Skin 

-4609.1 

523 

46 

0.0 

ext 

Upr.Skin 

-4924.3 

885 

46 

0.0 

ext 

Upr.Skin 

-4924.3 

885 

48 

-13.0 

ext 

Upr.Skin 

-4522.1 

47 
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Table I (Continued) 

Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

50 

0.0 

ext 

Upr.Skin 

-4465.8 

518 

51 

0.0 

ext 

Upr.Skin 

-4430.3 

877 

52 

45.0 

ext 

Upr.Skin 

-543.4 

877 

53 

90.0 

ext 

Upr.Skin 

1679.8 

877 

58 

0.0 

ext 

Upr.Skin 

-3050.7 

125 

62 

0.0 

ext 

Upr.Skin 

-3581.4 

872 

63 

0.0 

ext 

Upr.Skin 

-4809.6 

1328 

67 

0.0 

ext 

Upr.Skin 

-1910.8 

121 

65 

0.0 

ext 

Upr.Skin 

-4257.0 

864 

67 

0.0 

ext 

Upr.Skin 

-1910.8 

121 

69 

0.0 

int 

Int.Coast 

1201.7 

Q3098/99 

70 

45 “1 

int 

Int.Coast 

535.5 

Q3098/99 

71 i 

90 

int 

Int.Coast 

-681.0 

Q3098/99 

72 1 

0.0 

ext 

Upr.Skin 

-3473.1 

68 1 

74 

-13.0 

ext 

Upr.Skin 

-2434.0 

107 

75 

0.0 

ext 

Upr.Skin 

-6015.1 

117 

77 

-13.0 

ext 

Upr.Skin 

-1324.6 

221 

78 

0.0 

ext 

Upr.Skin 

-8430.2 1 

43 

79 

0.0 

ext 

Upr.Skin 

-7868.0 

25 

81 

0.0 

ext 

Upr.Skin 

-6647.8 

62 

84 

0.0 

ext 

Upr.Skin 

-5017.2 

63 

92 1 

0.0 

ext 

Upr.Skin 

-5082.1 

1234 

93 

0.0 

ext 

Upr.Skin 

-5485.3 

1292 

94 I 

0.0 | 

ext 

Upr.Skin 

-2647.8 

614 

95 

0.0 

ext 

Upr.Skin 

-2647.8 

614 

96 

0.0 

ext 

Upr.Skin 

-4978.8 

1359 

97 

0.0 

ext 

Upr.Skin 

-4978.8 

1359 1 
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Table I (Continued) 


Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

102 

0.0 

int 

Lwr.Skin 

2231.2 

1467 

123 

0.0 

int 

Lwr.Skin 

1957.8 

2244 

125 

0.0 

int 

Lwr.Skin 

2475.8 


128 

0.0 

int 

Lwr.Skin 

2616.2 

2276 

136 

0.0 

int 

Lwr.Skin 

2493.4 

1371 


0.0 

int 

Lwr.Skin 

3138.0 

1452 

1 144 

0.0 

ext 

Lwr.Skin 

3184.9 

1605 

■■m 


aft 


-280.4 

3007 



aft 


-139.7 

2955 

209 

-45.0 

aft 


-664.6 

2891 

302 

0.0 

aft 


-1102.3 

3324 

303 

45.0 

aft 


27.0 

3324 

304 

90.0 

aft 


352.9 

3324 

306 

0.0 

aft 


- 779.5 


307 

45.0 

aft 


136.2 


308 

90.0 

aft 


348.4 

3307 

310 

0.0 

aft 


-194.1 

3292 

311 

45.0 

aft 



3292 

312 

90.0 

aft 

Rear.Spar 

-383.5 

3292 

444 

45.0 

inb 

RIBA 

-380.3 

3703 

445 

90.0 

inb 

RIBA 

-114.1 

3651 

448 

0.0 

inb 

RIBA 

986.5 

3674 

454 

45.0 

inb 

RIB. 5 

530.0 

3918 

455 

90.0 

inb 


56.2 

3823 

464 

-45.0 

inb 


168.7 

3960 

465 

90.0 

inb 


50.6 

3961 
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Table I (Continued) 

Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

468 

0.0 

inb 

Rib.6 

-318.3 

4167 

474 

-45.0 

inb 

Rib.7 

-173.0 

4348 

475 

90.0 

inb 

Rib.7 

92.6 

4173 

478 

0.0 

inb 

Rib.7 

-195.7 

4250 

484 

-45.0 

inb 

Rib. 8 

-180.8 

4475 

485 

90.0 

inb 

Rib.8 

82.1 

4382 

488 

0.0 

inb 

Rib. 8 

549.8 

1 4507 

501 

0.0 

fwd 

Ext.Box 

138.6 

5104 

502 

-45.0 

fwd 

Ext.Box 

-2483.7 

5104 

503 

90.0 

fwd 

Ext.Box 

-316.5 

5104 

507 “ 1 

0.0 

int 

Ext.Box 

282 

Q4160 

511 

0.0 

ext 

Tr.Str.Up 

55.1 

4852 

513 

oxT 

ext 

Tr.Str.Up 

47.3 

4865 

602 

0.0 "I 

ext 

Upr.Skin 

-4670.7 

387 

604 1 

0.0 

ext 

Upr.Skin 

-4924.9 

421 

606 

0.0 1 

ext 

Upr.Skin 

-5083.4 

453 

610 

-13.0 

ext 

Upr.Skin 

-1618.8 

221 

611 

0.0 | 

ext 

Upr.Skin 

-3885.6 

49 

612 

0.0 

ext 

Upr.Skin 

-8596.0 

39 

613 | 

0.0 

ext 

Upr.Skin 

-2281.1 

120 

614 

0.0 

ext 

Upr.skin 

-2489.0 

69 

615 

0.0 t 

ext 

Upr.skin 

-4350.3 

444 

616 

0.0 

ext 

Upr.skin 

-4551.4 

449 

617 

0.0 

ext 

Upr.skin 

-4924.3 

885 
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Table II. Strain prediction from global model at ultimate load of 
166,000 lbs (including all gages on internal upper skin). 


Strain 
Gage No 


Surface I Location 



Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Micro 

Strain 


-4015.2 

-4488.6 

-4780.7 

-4780.7 

-4389.4 


-4389.4 


1607.8 


-4032.7 


400.9 

520.7 

708.9 
-4242.4 
-4326.3 

500.4 

279.8 


595.0 
-4986.1 
-2526.6 

1627.0 
-2639.8 
-3500.6 
-5589.5 
-5071.2 


-5176.1 


2582.9 


2803.8 


496.4 


-1096.2 


Node 
Number 
in FEM 
Model 


196 

194 


1359 


1359 


614 



474 


523 

4388 

4388 


4388 

885 


885 


493 

1328 


121 

25 

1292 

1234 


496 

550 


1550 


1550 
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Table II (Continued) 

Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 


Angle 


Surface 


Location 


Micro 

Strain 


107 

109 

111 

112 

113 

115 


0.0 

0.0 

0.0 

45.0 

90.0 
0.0 


ext 

ext 


Lwr.Skin 


Lwr.Skin 


ext 


Lwr.Skin 


ext 

ext 


Lwr.Skin 


Lwr.Skin 


ext 


Lwr.Skin 


2640.0 


3342.7 


2898.7 


975.5 

-1076.0 


3001.5 


116 


0.0 


ext 


Lwr.Skin 


2535.1 


118 


119 

120 
122 
124 


126 

129 

131 

132 

133 

137 

138 

139 
202 

203 

204 
206 

207 

208 
210 
211 
212 

213 

214 


0.0 


45.0 

90.0 

13.0 
13.0 


0.0 

0.0 

0.0 

45.0 

90.0 


13.0 
0.0 
0.0 
0.0 

-45.0 

90.0 
0.0 

-45.0 

90.0 

0.0 

-45.0 

90.0 

0.0 

0.0 


ext 


ext 


ext 


ext 


ext 


ext 


ext 


ext 


ext 


ext 


ext 

ext 


ext 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 

fwd 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 

Fwd.spar 


2296.7 

570.7 


-706.9 


794.6 


2099.3 

2354.4 


2627.3 


2620.1 


405.3 


-843.0 


3169.2 

2936.5 


3050.8 

-673.2 

-370.2 

-574.1 

-541,3 

-107.8 

98.9 

-340.4 

6.1 

43.0 

-3958.6 

-4064.6 


Node 
Number 
in FEM 
Model 

1667 

1558 

1574 

1574 

1574 

1446 

1691 

1784 

1784 

1784 

1385 



/U 

2291 

3007 

3007 

2007 

2955 

2955 

2955 

2891 

2891 

2891 

3172 

2671 


14 
















































































































Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 


215 


216 


217 


218 

301 

305 

309 

441 

442 

443 
449 

451 

452 


453 

461 

462 

463 
469 

471 

472 

473 

479 

480 


481 

482 

483 
489 
504 

504 

505 
505 


506 


Angle 


Surface I Location 


45.0 


0.0 

45.0 


90.0 

0.0 


0.0 

45.0 


90.0 


0.0 


-45.0 

90.0 


0.0 

0.0 


-45.0 


90.0 

0.0 


90.0 


0.0 

-45.0 


90.0 


0.0 

0.0 

0.0 

-45.0 

-45.0 


90.0 




RIBA 


RIBA 


RIBA 

RIBA 


RIB. 5 
RIB. 5 


RIB. 5 


Rib.6 


Rib.6 


Rib.6 

Rib.6 

Rib.7 


Rib.7 


Rib.7 


Rib.7 


Rib.7 


Ri 


Ri 


Rib.8 


Rib. 8 


Ext.Box 


Ext.Box 


Ext.Box 


Ext.Box 


Ext.Box 


Micro 

Strain 


-3585.7 

-3238.2 

-2217.4 

-1293.1 

66.3 


75.5 


1023.8 

267.7 
-422.5 
-260.4 

998.7 
-185.5 
532.0 


292.7 


133.6 

184.2 

-108.5 

-549.9 

-303.4 


-149.6 


211.8 


-299.5 


461.0 


0.0 


-2456.7 


0.0 


-319.2 


Node 
Number 
in FEM 
Model 


2815 


2836 


3324 


3307 

3292 

3703 

3703 

3703 

3674 

3918 

3918 

3918 


3960 

3960 

3960 

4167 

4348 


4348 


4348 


4250 


4334 


5104 


5104 


5104 


5104 


5104 


15 

























































































Table II (Continued) 

Strain prediction from global model at ultimate load of 166,000 lbs. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

506 

90.0 

aft 

Ext.Box 

0.0 

5104 

TO 

0.0 

int 

Tr.Str.Up 

66.3 

4865 

601 

0.0 

int 


-4472.7 

387 

609 

-13.0 

int 

ES3S3591 

-4525.4 

r 57 

141 

0.0 

ext 

Lwr.Skin 

3338.4 

1452 

144 1 

0.0 | 

ext 

Lwr.Skin 

3288.1 

1605 

146 1 

0.0 

ext 

Lwr.Skin 

3220.3 

1636 

148 

0.0 

ext 

Lwr.Skin ! 

3136.2 

1674 


0.0 

ext 

Lwr.Skin 

3038.0 

1714 

152 

0.0 1 

ext 

Lwr.Skin 

2926.0 

1753 

154 

0.0 

ext 

Lwr.Skin 

2796.6 


156 

0.0 

ext 

Lwr.Skin 

2623.0 

1836 


16 






















































Table III. Srain prediction from global 
all gages on external upper s\ 


Strain 
Gage No 

Angle 

Surface 

2 

0.0 

ext 

4 

0.0 

ext 

5 

0.0 

ext 

6 

45.0 

ext 

7 

90.0 

ext 

11 

0.0 

ext 

12 

0.0 

ext 

15 

0.0 

ext 

16 

0.0 

ext 

17 

0.0 

ext 

18 

0.0 

ext 

20 

0.0 

ext 

22 

0.0 

ext 

25 

0.0 

ext 

26 

0.0 

ext 

27 

45.0 

ext 

28 

90.0 

ext 

32 

0.0 

ext 

37 

0.0 

ext 

39 

0.0 

ext 

46 

0.0 

ext 

46 

0.0 

ext 

48 

-13.0 

ext 

50 

0.0 

ext 

51 

0.0 

ext 

52 

45.0 

ext 

53 

ffff 1 ' 1 ' NffP 

ext 

58 


ext 


C C 


at 50% of ultimate load (including 



Location 






Micro 

Strain 

Node 
Number 
in FEM 
Model 

-1805.0 

196 

-1915.2 

194 

-2428.3 

65 

-827.9 

65 

915.4 

65 

-2430.7 

1359 

-2430.7 

1359 

-2268.3 

614 

-2268.3 

614 

-2239.3 

71 

-2339.6 

414 

-2266.2 

437 

-2105.2 

371 

-2113.0 

101 

-2195.9 

403 

-397.5 

403 

703.7 

403 

-2019.1 

401 

-2299.2 

474 

-2258.8 

523 

-2416.9 

885 

-2416.9 

885 

-2242.0 

47 

-2342.2 

518 

-2305.2 

877 

-376.4 

877 

738.7 

877 

-1934.7 

125 






































































Table III (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

62 

0.0 

ext 


IHEZ&S9 

872 

63 

0.0 

ext 

ilkiigawsM 

■na 

1328 

65 

0.0 

ext 


-1869.7 

864 

67 


ext 


-1518.6 

121 

69 


int 

Int.Coast 

596.8 

03098/99 

70 J 

45 

int 

Int.Coast 

430.6 

03098/99 

71 J 


int 

Int.Coast 

-257.5 

Q3098/99 

72 

0.0 

ext 



68 

74 ! 

-13.0 

ext 

WB3M3M 


107 

75 

0.0 

ext 

KBWSBBB 


117 

77 1 

-13.0 

ext 


-810.5 

221 

78 

0.0 

ext 

Ufflyfi/tH 

1WiCT| 

43 

79 

0.0 

ext 


ma 

25 

81 

0.0 

ext 

Biii'Ti him 

-2889.8 

62 

84 ! 

0.0 1 

ext 


-2308.8 

63 

92 

0.0 1 

ext 


-2545.0 

1234 

93 

0.0 

ext 


-2751.8 

1292 

94 ] 

0.0 1 

ext 



614 

95 

o.o I 

ext 


■Kfetalsl 

614 

96 

0.0 

ext 

■ssssn 

-2430.7 

1359 

102 1 

0.0 

int 

Lwr.Skin 

1113.9 

1467 

123 

0.0 } 

int 

Lwr.Skin 

982.1 

2244 


18 



















































































Table III (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

125 

0.0 

int 

Lwr.Skin 

1238.6 

2260 

128 

0.0 

int 

Lwr.Skin 

1302.9 

2276 

136 

0.0 

int 

Lwr.Skin 

1253.9 

1371 

140 

0.0 

int 

Lwr.Skin 

1554.2 

1452 

144 

0.0 

ext 

Lwr.Skin 

1586.5 

1605 

201 

-45.0 

aft 


-117.9 

3007 

205 

-45.0 

aft 


-63.4 

2955 

209 

-45.0 

aft 


2891.0 

302 

303 

45.0 

aft 


21.2 

3324 

304 

90.0 

aft 


176.3 

3324 

306 

0.0 

aft 


-398.8 

3307 

307 

45.0 

aft 


60.0 

3307 

308 

90.0 

aft 


170.3 

3307 

310 

0.0 

aft 


-94.9 

3292 

311 

45.0 

aft 


518.7 

3292 

312 

90.0 

aft 


-198.3 

3292 

444 

45.0 

inb 

RIBA 

-189.3 

3703 

445 

90.0 

inb 

RIBA 

15.0 

3651 

448 

0.0 

inb 

RIBA 

492.7 

3674 

454 

45.0 

inb 

RIB.5 

289.8 

3918 

455 

90.0 

inb 

RIB. 5 

43.6 

3823 

464 

-45.0 

inb 

Rib.6 

111 A 

3960 

465 

90.0 

inb 

Rib.6 

57.2 

3961 

468 

0.0 

inb 

Rib.6 

-174.2 

4167 

474 

-45.0 

inb 

Rib.7 

-39.2 

4348 

475 

90.0 

inb 

Rib.7 

65.9 

4173 

478 

0.0 

inb 

! Rib.7 

-101.9 

4250 

484 

-45.0 

inb 

Rib.8 

-105.1 

4475 

485 

90.0 

inb 

Rib.8 

61.7 

4382 j 

488 

0.0 

inb 

Rib.8 

328.2 

4507 | 

501 

0.0 

fwd 

Ext.Box 

79.5 

5104 


19 















































































































































































Table III (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

502 

-45.0 

fwd 

Ext.Box 

-1240.2 

5104 

503 

90.0 

fwd 

Ext. Box 

-157.4 

5104 

507 

0.0 

Int 

Ext.Box 

161.6 

Q4160 

511 

0.0 

ext 

Tr.Str.Up 

27.0 

4852 

513 

0.0 

ext 

Tr.Str.Up 

23.2 

1 4865 

602 

0.0 

ext 


-2317.1 

387 

604 

0.0 

ext 


-2322.1 

421 

606 

0.0 

ext 


-2191.3 

453 

610 

-13.0 

ext 

mimvm 

-810.5 

221 

611 

0.0 

ext 


-1995.0 

49 

612 

0.0 

ext 


-4073.9 

39 

613 

0.0 

ext 


-1626.0 

120 

614 

0.0 

ext 

MBBW 

-1795.9 

69 

615 

0.0 

ext 



444 

616 

0.0 

ext 


wiw 

449 

617 

0.0 

ext 


-2416.8 

885 


20 





























































Table IV. Strain prediction from global model at 50% of ultimate load 
(including all gages on internal upper skin). 


Strain 
Gage No 


Angle 





0.0 

oo 

0.0 


45.0 


90.0 


0.0 




Surface Location 








Upr.Skin 

Upr.Skin 

Upr.Skin 

Upr.Skin 

Upr.Skin 


U 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Micro 

Strain 


-2001,4 

-2206.4 

-2358.6 

-2358.6 

-2172.0 

-2172.0 


-2119.2 


-2014.1 

137.0 

154.1 

259.7 
-2134.3 
-2185.8 

245.7 

109.7 

298.1 
-2535.9 


-2535.9 

890.6 


-1761.1 

-1613.7 

-2970.6 


-2547.4 


-2570.4 


1291.5 


1402.1 


247.8 


Node 
Number 
in FEM 
Model 


196 

194 

1359 

1359 

614 


614 


399 


401 


4336 


474 

523 

4388 

4388 

4388 


493 


1328 

121 

25 

1292 


1234 


1496 


1550 


1550 


21 



















































































Table IV (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 



Angle 


90.0 


0.0 


0.0 


0.0 


45.0 


90.0 


0.0 


0.0 


0.0 


45.0 


90.0 


13.0 


13.0 


0.0 


0.0 


0.0 


45.0 


90.0 


13.0 


0.0 


0.0 


0.0 


-45.0 


90.0 

0.0 


-45.0 


90.0 


0.0 


Surface I Location 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Lwr.Skin 


Micro 

Strain 


-540.1 


1450.6 


490.9 


-522.7 


1486.0 


290.1 


-361.1 


312.2 


694.3 


1181.3 


315.3 

323.5 


217.2 


1469.8 


Node 
Number 
in FEM 
Model 


1550 


1667 


1558 


1574 


1574 


1574 


1446 


1691 


1784 


1784 


1784 


1385 


1374 


2273 


2276 



22 








































































































Table IV (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

211 

-45.0 

fwd 


17.9 

2891 

212 

90.0 

fwd 


15.4 

2891 

213 

0.0 

fwd 



3172 

214 

0.0 

fwd 

H33SES91 


2671 

215 

0.0 

fwd 



2719 

216 

0.0 

fwd 



2767 

217 

0.0 

fwd 


-1087.3 

2815 

218 

0.0 

fwd 


-618.9 

2836 

301 


fwd 


34.1 

3324 

305 


fwd 


54.8 

3307 

309 

45.0 

fwd 


521.4 

3292 

441 

0.0 

out 

RIB .4 

125.3 

3703 

442 

45.0 

out 

RIB .4 

-196.2 

3703 

443 

90.0 

out 

RIBA 

-32.2 

3703 

449 

0.0 

out 

RIBA 

515.8 

3674 

451 


out 

RIB. 5 

-89.3 

3918 

452 


out 

RIB. 5 

294.0 

3918 

453 

90.0 

out 

RIB. 5 

170.6 

3918 

461 

0.0 

out 

Rib.6 

84.5 

3960 

462 

-45.0 

out 

Rib.6 

111.4 

3960 

463 

90.0 

out 

Rib.6 

-39.4 

I 3960 

469 

0.0 

out 

Rib.6 

-300.4 

4167 

471 

0.0 

out 

Rib.7 

- 77.0 

4348 

472 

-45.0 

out 

Rib.7 

-34.7 

4348 

473 

90.0 

out 

Rib.7 

106.8 

4348 

479 

0.0 

out 

Rib.7 

-161.3 

4250 

480 

90.0 

out 

Rib.7 


4334 

481 

0.0 

out 

Rib.8 

- 57.2 

4475 






















































































Table IV (continued) 


Strain prediction from global model at 50% of ultimate load. 


Strain 
Gage No 

Angle 

Surface 

Location 

Micro 

Strain 

Node 
Number 
in FEM 
Model 

482 

-45.0 

out 

Rib.8 

-95.9 

4475 

483 

90.0 

out 

Rib.8 

154.2 

4475 

489 

0.0 

out 

Rib.8 

569.8 

4507 


0.0 

aft 

Ext.Box 

86.1 



0.0 

aft 

Ext.Box 

0.0 

■Em 

505 

-45.0 

aft 

Ext.Box 

-1233.4 

■ i 

505 

-45.0 

aft 

Ext.Box 

0.0 

■E B 

506 

90.0 

aft 

Ext.Box 

-158.9 

| 

506 

90.0 

aft 

Ext.Box 

0.0 

| 

510 

0.0 

int 

Tr.Str.Up 

32.6 


601 1 

0.0 

int 

BflSSiSSi 

-2222.2 

387 

609 

-13.0 

int 

ESSE 

-2495.9 

57 

141 

0.0 

ext 

Lwr.Skin 

1654.5 

1452 

144 

0.0 

ext 

Lwr.Skin 

1648.8 

1605 

146 

0.0 

ext 

Lwr.Skin 

1615.9 

1636 

148 

0.0 

ext 

Lwr.Skin 

1573.7 

1674 

150 

0.0 

ext 

Lwr.Skin 

1524.6 

1714 

152 

0.0 

ext 

Lwr.Skin 


1753 

154 

0.0 | 

ext 

Lwr.Skin 


1800 1 

156 

0.0 

ext 

Lwr.Skin 

1315.9 

1836 


24 


























































































Table V Strain estimation from local model at ultimate load 
of 166,000 lbs. 
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Table VI Strain estimation from local model at 83,000 lbs (50% 
Loading) 


Strain 
Gage No 

Angle 

Surface 

Micro 

Strain 

Node 

22 

0 

Ext 

-2533.00 

18 

23 

0 

IntINB 

-22.50 

E2000 

24 

0 

IntOUT 

-62.50 

E2000 I 

33 

0 

Int 

1 457.00 

E2415 

34 

0 

Int 

-442.50 

E2415 

35 

0 

Int 

-480.00 

E2415 

43 

0 

Int 

-2454.75 

1184 

44 

r o 

Int 

-2784.50 

1203 

45 

r o 

Int 

-2585.00 

1214 

46 

0 

Ext 

-3292.00 

1212 

48 

-13 

Ext 

-2807.77 

239 

49 

0 

INT 

-1231.50 

506 

50 

0 

Ext 

-2275.25 

504 

51 

0 

Ext 

-2277.25 

402 

52 

45 

Ext 

-1212.63 

402 

53 

90 

Ext 

694.25 

402 

54 

0 

INT 

-1436.50 

506 

55 

0 

Int 

-2670.00 

2643 

56 

0 

Int 

-3227.00 

2613 

57 

0 

Int 

-2928.00 

751 

58 

0 

Ext 

-2248.50 

498 

59 

0 

Int 

-2150.25 

1886 

60 

0 

Int 

-2553.50 

676 

61 

0 

Int 

-2434.25 

1891 

62 

0 

Ext 

-2276.00 

679 

63 

0 

Ext 

-2271.25 

1512 

64 

0 

Int 

-1999.50 

1512 

65 

0 

Ext 

-2148.50 

1296 

66 

0 

IntINB 

-867.50 

E2130 

r 607 

0 

IntINB 

-332.00 

E2070 

608 

0 

Ext 

-2246.25 

468 

609 

-13 

Int 

-1665.24 

242 

617 

0 

Ext 

-2867.50 

2148 


26 
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Figure 1 . Photograph of McDonnell Douglas wing box test article. 
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Figure 2. Dimensions of McDonnell Douglas wing box test article. 


Outboard 
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Figure 3. Upper skin panel construction. 







Thickness (Number of Laminate Stacks) 
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Figure 4. Stub box upper cover thickness distribution. 
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Figure 6. Finite element model of McDonnell Douglas composite wing stub box. 
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Figure 7. Interior of the stub box finite element model. 
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Figure 8. Local finite element model of wing stub box. 
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Figure 10. Details of stiffener runout region. 
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Figure 1 1 . Buckling mode of composite wing box. 
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Figure 14. Deformed composite box model due to 166,000 lb. tip load. 


W along Line A-A 
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Figure 15. Top cover out-of-plane displacements due to 166,000 lb. tip load. 



c 
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Load Factor 

Figure 16. Lower skin out-of-plane displacements at instrumented locations due to 166,000 lb. tip load. 



43 


Figure 1 7. Predicted longitudinal strain contours at skin-side surface of upper cover 
due to 166,000 lb. tip load. 



inVin. 
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Figure 18. Predicted longitudinal strain contours at middle surface of upper cover 
due to 166,000 lb. tip load. 
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Figure 19. Predicted longitudinal strain contours at stiffener-side surface of upper cover 
due to 166,000 lb. tip load. 
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Figure 20. Predicted contours of strain circumferential to hole at skin-side surface of upper cover 
due to 166,000 lb. tip load. 





-0.006770 
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Figure 21 . Predicted contours of strain circumferential to hole at mid-surface of upper cover 
due to 166,000 lb. tip load. 
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Figure 22. Predicted contours of strain circumferential to hole at stiffener-side surface 
of upper cover due to 166,000 lb. tip load. 
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Figure 23. Upper skin strain gage locations. 
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Figure 24. Lower skin strain gage locations. 
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Figure 26. Predicted strains at locations of gage 78. 
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Figure 27. Predicted strains at location of gage 79. 




54 


Figure 28. 
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Figure 29. Predicted strains at locations of gages 63 and 64. 
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Strains (microstrain) 

Figure 30. Predicted strains at the locations of gages 67 and 68. 
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Figure 31 . Predicted deformed shape of local model due to 166,000 lb. tip load. 


-0.006120 





58 


Figure 32. Predicted longitudinal strain contours at skin-side surface 
of upper cover due to 166,000 lb. tip load. 
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Figure 33. Predicted longitudinal strain contours at mid-surface of upper cover due to 
166,000 lb. tip load. 
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Figure 34. Predicted longitudinal strain contours at stiffener-side surface of upper cover 
due to 166,000 lb. tip load* 
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Figure 35. Predicted longitudinal strains for gage at runout 1. 
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Figure 36. Predicted longitudinal strains for gage at runout 2. 
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Figure 37. Predicted longitudinal strains for gage at runout 3. 







Figure 39. Predicted lateral strains on mid-surface of the ribs in local model 
due to 166,000 lb. tip load. 
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Figure 40. Predicted lateral strains on outboard surface of ribs in local model -0.004010 

due to 166,000 lb. tip load. 
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Figure 41 . Predicted normal stress resultant (separation load) between stiffener 1 and skin in upper cover. 
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X Location 

Figure 42. Predicted normal stress resultant (separation load) between stiffener 2 and skin of upper 
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Figure 43. Predicted normal stress resultant (separation load) between stiffener 3 and the skin of the upper cover. 
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Figure 44. Predicted normal stress resultant (separation load) between stiffener 4 and the skin of the 
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